Symmetry, dimension and the distribution of the conductance at the mobility edge 
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The probability distribution of the conductance at the 
mobility edge, p c (g), in different universality classes and di- 
mensions is investigated numerically for a variety of random 
systems. It is shown that p c (g) is universal for systems of 
given symmetry, dimensionality, and boundary conditions. 
An analytical form of p c (g) for small values of g is discussed 
and agreement with numerical data is observed. For g > 1, 
In Pc(g) is proportional to (g — 1) rather than (g — l) 2 . 
PACS numbers: 71.30.+h, 71.55.Jv 
Disordered systems may show a transition frcffla,metal- 
lic states to insulating ones at a mobility edge,Em which 
seperates the two regions. The probability distribution 
p(g) of the conductance g at the critical point of dis- 
ordered systems undergoing a localizatiea-delocalization 
transition is still under investigationHliiJ Previous stud- 
ies have shown that p{g) depends on the symmetry of 
the system,cl its dimensionality ,0 the boundary .conditions 
perpendicular to the direction of transport ,cnlj and the 
amount of anisotropy.B'tL 2 ! Yet, a complete theoijY-fixplain- 
ing the form of the distribution is still missingoU Know- 
ing that the conductance distributions are normal and 
log-normal in the extended, metallic and the localized, 
insulating regimes respectively, and taking into account 
the continuous nature of the Anderson localization- 
delocalization transition, it seems reasonable to try to 
combine the two forms. 

In this letter we have calculated the probability distri- 
bution of the conductance for a variety of systems of dif- 
ferent dimension and symmetry, give an approximate ex- 
pression for p c {g)i compare our numerical results to some 
analytical approximations, and present a way of explain- 
ing the differences between theory and numerical results. 
We find that it is necessary to examine the distributions 
of the smallest Lyapunov exponents and the relationship 
between their respective mean values. We also show that 
p c {g) is independent of the particular point chosen on the 
critical surface in parameter space, consistent with sim- 
ilar findings on v aryin g the distribution function of the 
disorder potential.HEJ 

We have calculated the conductance distributions at 
the mobility edge of a three-dimensional (3D) system 
with orthogonal symmetry, a two-dimensional (2D) sys- 
tem with symplectic symmetry, and of several 2D systems 
with unitary symmetry. All these systems possesss a mo- 
bility edge and are modelled after the Anderson tight- 
binding Hamiltonian 
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lattice. The variables r, r' take on values of 1 or —1 for 
symplectic systems with spin-orbit interactions, where 
the hopping integrals V n y thus become 2x2 matrices; 
otherwise they .are scalars and the spin "variables" have 
only one value .113 The site energies e n are always inde- 
pendent of r. 

In Fig. [j] we show three unitary systems with peri- 
odic boundary conditions. A magnetic field perpendicu- 
lar to the direction of transport facilitates the existence 
of critical states at the center of each Landau subband. 
We investigate the dependence of p c (g) on the disorder 
strength W. The flux per unit area a has been kept con- 
stant at a = 1/8. For weak disorder, considerable finite 
size effects have to be eliminated. Even for the system 
shown with 192 x 192 lattice sites (dashed line in Fig. 
|l|) the distribution still has not completely converged to 
the form obtained for the two cases of stronger disorder, 
where the system size is only 64 x 64 lattice sites (solid 
lines in Fig. [j]). Anisotropic systems can be rescaled to 
the same distribution.t 2 ] Table Q contains the averages and 
standard deviations of the relevant variables for these en- 
sembles as well as for those we use in later parts of this 
paper. 

We will discuss the transmission properties of a system 
in terms of its "extensive Lyapunov exponents" where 
e Zi are the eigenvalues of T^T and T is the transfer matrix 
in the channel representation, p3jh.cn, we have for the 
conductance g (in units of e 2 /h)£3 
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where N is the number of open channels. The distribu- 
tion of the conductance should therefore be discussed in 
connection with that of the Lyapunov exponents. The 
distribution of the smallest positive Lyapunov exponent 
z\ can be approximated by a Wigner distribution with 
= 1, independently of the actual universality class of 
the system:BEZI 
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where n, n' are nearest neighbour sites in the 2D or 3D 



where (.) denotes the ensemble average. This approxima- 
tion works reasonably well, if (zi) is small enough, which 
is true in two and three dimensions, but not e.g. in four. 

Approximating g i=s g\ A = cosh~ 2 (zi/2), we can rewrite 
this distribution in terms of ln(g) as 

p(ln g) «p(ln 51) 
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evaluated at \n.(g) = — 2 lncosh(zi/2). This obviously 
neglects contributions to the conductance from higher 
channels and therefore overestimates the distribution in 
the range hi(g) < 0. Note, that, because cosh 2 (zi/2) > 1 
for all zi, ln(gi) < 0. One finds that p(\n gi) is already 
in reasonable agreement with p(hi g) indicating that the 
higher channels' contributions are small, though not en- 
tirely negligible. Therefore p(ln g\) can be used as a 
starting point for discussion of the correct distribution 
of the conductance in the range g < 1. Fig. |2| shows the 
numerical results for 10, 000 cubic systems of orthogo- 
nal symmetry with periodic boundary conditions. It can 
be seen clearly from Fig. || that both p(ln g{) and Eq. [| 
are in very good agreement with the detailed numerical 
results. Also shown is the distribution p(ln g2), where 
.92 = gi + cosh~ 2 (z2/2), which agrees already very well 
with the distribution of the total conductance. Squares 
of symplectic symmetry behave similarly. Also, systems 
with hard wall boundary conditions show the same qual- 
itative behaviour in both 3D orthogonal systems and 2D 
symplectic ones. A summary of the averages and vari- 
ances of Z\ and zi can be found in Table |J. 

Using a different, more elaborate approach, Muttalib 
and WolffeM derived for quasi-one-dimensional, weakly 
disordered systems a formula for the critical probability 
distribution over the whole range of g, including g > 1. 
It can be written as 

{1 Vzisinhzi . < 1 

where the formula for the range g < 1 again needs to be 
evaluated at ln(g) = — 21ncosh(zi/2). The parameter T 
can be used to fit this function to the numerical results, 
(a is a function of T.) Taking T = n/ (zi) and noting 
that sinhfzi) « z\ for small Zi, the similarity of Eq. |5| 
and Eq. ^ is apparent. This suggests replacing z\ in the 
prefactor of Eq. || with yj z\ sinh(zi). Preliminary re- 
sults show that this actually results in better agreement 
with data even for somewhat higher values of (zi). It 
should be noted though, that in T instead of the average 
value of zi, one should use the most probable one, which 
is smaller than the average value by a factor of about 
0.8 in the case of a Wigner distribution. Despite their 
deriving^ a distribution for the whole range of g, their 
formula still overestimates slightly the weight of the range 
g < 1 in 3D systems. However, in the 2D symplectic 
case, they slightly underestimate this weight. A change 
in the fitting parameter T does not remedy this discrep- 
ancy in a satisfactory manner. Fig. || shows numerical 
results together with a fit according to Eq. [jj. The first 
panel shows the distributions for 3D orthogonal systems 



with 10 x 10 x 10 lattice sites. In a log-linear plot one 
can see that Eq. || increasingly overestimates p(ln g) for 
ln(<jf) — > —00. For the 2D symplectic case shown in the 
second panel, a fit for In (5) close to results in a very 
strong underestimation far from ln(<?) = 0. The fit pre- 
sented for both kinds of boundary conditions still gives an 
overall undexstimation of the portion of the conductance 
distribution!^ with g < 1. 

To understand the qualitatively different behaviours of 
this theoretical approach, one has toiook at the averages 
of the higher Lyapunov exponents £3 In the quasi-one- 
dimensional, weakly disordered case for which Eq. ^ was 
derived, one has (22} = 2 • (zi), indepen4wit|--of dimen- 
sion, symmetry or boundary conditions.E^nlZHlj For the 
3D orthogonal ensemble, one finds at the critical point 
that (zi^H oc i, and thus (22) is significantly smaller than 
2 • (^i)J— 1 so that the contribution of the second chan- 
nel is higher than expected from Eq. whereas for the 
2D symplectic case, (22) > (zi), so that the second chan- 
nel's contribution is smaller than expected. Compare the 
values in Table O, which support these arguments. 

We also looked at the conductance distribution in the 
range g > 1. In order to have a sizeable ensemble for that 
range, we took half a million samples for 2D symplectic 
systems of 40 x 40 lattice sites with both periodic and 
hard wall boundary conditions as well as for 3D orthog- 
onal systems of 10 x 10 x 10 lattice sites with periodic 
boundary conditions. For cubes with hard wall bound- 
ary conditions we even took ten million samples. About 
20% of the symplectic samples, 6% of the 3D samples 
with periodic boundary conditions, and 3% of the 3D 
samples with hard wall boundary conditions turn out to 
have a conductance bigger than 1. For the latter ensem- 
ble, only about 470 out of the ten million samples have 
a conductance g > 2 and only one sample can be found 
with g > 3. We find that in the range g > 1, In p c (g) 
is at most linear in {g — 1), as can be seen from Fig. 
f§ This is in disagreement with the theory presented by 
Muttalib and W61fle,Eil which predicts a quadratic depen- 
dence with a logarithmic correction, and which therefore 
expects a positive first derivative of h\p c (g) in g. Fi- 
nally, Fig. ^ shows that the first derivative of p c (g) is 
discontinuoust3 at g = 1. We suppose that this non- 
analytical behaviour was not taken into account by the 
analysis of Muttalib and Wolfle.EJ 

In conclusion, we have shown that the critical distribu- 
tion of the conductance in disordered systems is univer- 
sal for systems of given dimensionality, universality class, 
and boundary conditions. We show further that for sys- 
tems of quite different types, the total conductance is 
distributed only slightly differently from the distribution 
of the first channel, and give arguments for the quality 
of corrections depending on the statistics of the second 
channel. We present a formula for p c (ln g) which agrees 
reasonably well with the numerical results in the range 
g < 1. Finally, we found non-analycity of p c (g) at g ss 1 
and estimated an exponent of roughly 1 in In p c (g) as a 
function of g — 1 rather than the predicted exponent of 
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FIG. 1. The conductance distribtuions for three different 
critical two-dimensional systems with a magnetic field per- 
pendicular to the plane. 
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FIG. 2. The distribution of the total conductance g of 
cubes with 10 x 10 x 10 lattice sites and periodic boundary 
conditions, together with the distributions for the contribu- 
tions from the first (g\) and the first two (32) channels of the 
same ensemble. The thin solid line is the result of Eq. ^ with 
(zi) = 2.825. 
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FIG. 3. The conductance distributions for a) 

three-dimensional systems of orthogonal symmetry with 
10 x 10 x 10 lattice sites and b) two-dimensional systems of 
symplectic symmetry with 40 x 40 lattice sites (thick lines). 
The thin lines are fits to the data according to Eq. ^ with 
T = tt/(2.2) 2 , T = tt/(2.55) 2 , and F = tt/(1.4) 2 for the 3d 
periodic, 3d hard wall, and 2d cases respectively. 
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FIG. 4. The distribution for an ensemble of 500, 000 cubic 
systems (thick line) shows a behaviour In p — const. + (g — l) a 
with q w 1 (thin line) in the region g > 1 and a discontinuity 
in its first derivative at g = 1. The distribution for an ensem- 
ble of 500, 000 square systems behaves similarly. Both cases 
shown use periodic boundary conditions. 

TABLE I. The averages and variances of the conductance 
and its logarithm for the ensembles we used in this work. O, 
U, S: orthogonal, unitary, symplectic; p, h: periodic, hard wall 
boundary conditions. Unitary systems use periodic boundary 
conditions. ^Vgtat- number of samples. 
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2d U, W = 4 


10,000 


0.445 


0.082 


-1.120 


0.842 


2d U, W = 2 


10,000 


0.428 


0.079 


-1.172 


0.887 


2d U, W = 0.5 


10,000 


0.393 


0.078 


-1.306 


1.027 


2d S, p 


500,000 


0.749 


0.088 


-0.401 


0.283 


2d S, h 


500,000 


0.691 


0.108 


-0.531 


0.418 


3d O, p 


500,000 


0.391 


0.108 


-1.418 


1.282 


3d O, h 


10,000,000 


0.284 


0.087 


-1.929 


1.762 



TABLE II. The averages and variances of the two smallest 
Lyapunov exponents. = 10, 000. 
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2d S, p 


1.424 


0.621 


3.987 


0.924 


2d S, h 


1.635 


0.811 


4.065 


1.186 


3d O, p 


2.825 


1.918 


4.965 


1.829 


3d O, h 


3.411 


2.475 


5.518 


2.132 
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